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In this Letter, we study the effect of topological zero modes on Entanglement Hamiltonians and 
entropy of free chiral fermions in (l+l)d. We show how Riemann-Hilbert solutions combined with 
finite rank perturbation theory allow us to obtain exact expressions for Entanglement Hamiltonians. 

In the absence of the zero mode, the resulting Entanglement Hamiltonians consists of local and bi¬ 
local terms. In the periodic sector, the presence of a zero mode leads to an additional non-local 
contribution to the entanglement Hamiltonian. We derive an exact expression for this term and for 
the resulting change in the entanglement entropy. 


Entanglement Hamiltonians (EH) are the next object 
to explore in a series of advances in our understanding 
the quantum structure of many-body states in field the¬ 
ory and condensed matter m- Indeed, much work has 
been devoted to understanding entanglement entropy, in 
a variety of systems, and, subsequently, more detailed in¬ 
formation about the spectrum of reduced density matri¬ 
ces is being investigated as well. Entanglement Hamilto¬ 
nians go one step further in that they contain information 
about the entanglement spectrum, as well as about the 
associated eigenvectors, and the possibility to understand 
the reduced state as a thermal state. Perhaps the most 
striking of the recent results about Entanglement Hamil¬ 
tonians is the realization that for a spherical entangling 
region in a conformal field theory (CFT), Entanglement 
Hamiltonians have a local form, that may be interpreted 
as the original hamiltonian energy density with a spa¬ 
tially dependent temperature EHU. This result is un¬ 
usual in its simplicity, as, unfortunately, the computation 
of EHs is, in general, much more involved than entropy 
and spectrum and only a few results are available. Thus 
it is a grand challenge to find additional solvable cases. 

In this Letter, we present a new method for comput¬ 
ing Entanglement Hamiltonians of free fermions in the 
presence of zero modes. Fermionic modes localized on 
topological defects often reflect the topological nature 
of the defect through unusual properties such as charge 
fractionalization [T5] and non-abelian braiding m and 
appear in a variety of systems from polyacetylene [14] 
to defects in topological insulators and superconductors 
USHHi- The modes are usually intimately tied to the 
appearance of ground state degeneracies of a topological 
nature. Here, we are able to present new results for the 
EHs for chiral fermions, and study in detail the effects of 
the boundary conditions (BC), periodic/anti-periodic for 
Majorana and generic for Dirac fermions, and of the zero 
modes (present in the case of periodic boundary condi¬ 
tions) on entanglement. The edge theory of the p + ip 
superconductor provides an explicit realization of such a 
scenario for Majorana fermions and serves as a concrete 
physical model for our calculation (see, e.g. [H]). Mean¬ 
while, our calculation in the Dirac sector is relevant to 


the edge theory of the Quantum Hall state. Wherever 
available we make contact with previous results in the 
literature. 

Our approach, as we will soon review, starts with 
the relation between the EH and the resolvent associ¬ 
ated with an appropriately projected (onto the entan¬ 
gling region) free fermion Green’s function, and adding 
the contribution of zero modes with an essentially exact 
re-summation. For concreteness we concentrate on chiral 
Majorana and Dirac fermions in (l + l)d, with the spatial 
direction taken to be a finite radius circle. As we find be¬ 
low, in this case, the computation of the resolvent can be 
recast into a Riemann-Hilbert problem (RHP, see below), 
generalizing the approach of Casini and Huerta pH who 
computed the EH for free chiral fermions, albeit without 
zero modes. We combine the RHP approach with ex¬ 
actly summable perturbation theory, to compute the EH 
for chiral Majorana and Dirac fermions on a finite circle. 
As our main result we find an exact analytic form for the 


entanglement Hamiltonian Eq. (16), with a contribution 


from the topological zero given separately in Eq. (22). 


The reduced density matrix py on a spatial region V is 
defined to reproduce expectation values of all operators 
Oy inside V via the relation Tr v (pvOv) = (Oy). The 
EH Hy is an effective Hamiltonian inside V given by 


Pv = 




Zy 


(1) 


For free fermions and in the absence of a ground state 
degeneracy (i.e. in the absence of zero modes), Wick’s 
Theorem implies that the EH is a quadratic operator 
whose kernel (the single particle EH) Hy is uniquely de¬ 
termined by the equal-time Green’s function 


G(x,y) = (iff(x)¥(y)) (2) 


via the relation Hy = — log ((PyGPy ) -1 — 1) [15] , Here 
PyGPy is the Green’s function restricted to the region V 
by the projectors Py. As we demonstrate in the supple¬ 
mentary materials, a similar relation holds for Majorana 
fermions ip, given (ip(x)ip(y)}, allowing us to treat the 
problem on a parallel footing. 
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The relation between Hy and G may be expressed in 
integral form as: 

/>00 

H v = -J i (L(0)+L(-0))dl3, (3) 

where L(/3) is the resolvent pQ 

L(/3) = (P V GP V -1/2+ p)-\ (4) 

Thus the derivation of the EH reduces to a computation 
of a projected Green’s function resolvent. It turns out 
that finding the resolvent L((3) can be recast as an RHP. 

This connection between the resolvent and the RHP 
arises in the following way: for free fermions (and in 
the absence of zero modes), the Green’s function G(x,y) 
computed on the ground state, im-restricted to some en¬ 
tangling region V, acts as a single particle operator that 
projects onto positive energy modes. For chiral fermions 
in (1 + 1 )d, the projection onto positive modes can be 
interpreted as a projection onto analytic functions on 
the upper/lower half of the complex plane. For such a 
G(x,y), the computation of the resolvent Q (associated 
with the projection of G onto V ) can be mapped onto a 
RHP (see, e.g. [2D]), which, in the simple case treated 
here is a scalar RHP and so exactly solvable. Succinctly 
put, a RHP is a jump problem along some smooth con¬ 
tour C in the complex plane, for a piecewise analytic ma¬ 
trix function. Solving it amounts to finding the matrix 
function X (z) which is analytic on C \C, and which has 
boundary values on both sides of C, X±(z) subject to a 
jump condition XZ 1 (s)2f + (s) = v(s) for s £ C, and with 
a given matrix v(s). The jump function for our case is a 
scalar determined by the projector onto V, Py. 

Our next idea is that once we have access to the resol¬ 
vent Q (either through the RHP or any other method), 
we can use it to explore various deformations of the sys¬ 
tem that yield simple perturbations of the Green’s func¬ 
tion Q. In particular, topologically nontrivial back¬ 
ground for the fermion will result in the appearance of 
zero modes Emm- As we show below, these manifest 
themselves as finite rank perturbations of the correla¬ 
tion matrices. Given the resolvent for the unperturbed 
Green’s function, we solve for the perturbed resolvent 
by summing up exactly the resulting perturbative series. 
Here we carry out this recipe for the simplest scenario of 
a chiral Dirac/Majorana fermion on a spatial circle with 
general boundary conditions. 

Starting with anti-periodic (Neveu-Schwarz, NS) 
fermions on a circle with no gauge field, inserting a n flux 
then takes the anti-periodic sector to the periodic (Ra- 
mond, R) sector, which has a zero mode whose effect on 
the EH we derive explicitly. Note that the appearance of 
ground state degeneracy for period boundary conditions 
has been studied extensively. In particular, Affleck and 
Ludwig [25], introduced the notion of boundary entropy, 
to understand the change in degeneracy due to impuri¬ 
ties. It was also noted in :24J that the boundary entropy 



FIG. 1: The entanglement region V = U (ajfij) for fermions on 
a ring of radius R considered in the paper, and an associated 
Riemann-Hilbert problem. 


can be manifested in a sub-leading contribution to the 
EE of a finite interval. While these calculations are usu¬ 
ally carried out by looking at the entire system, or at 
EE, our EH calculation allows us to study the effect of 
degeneracy from the point of view of the local density 
matrix. In Eq. ( |24| ), we find the exact change in the en¬ 
tanglement entropy (EE) of free chiral fermions between 
periodic and anti-periodic BC, for an arbitrary subset of 
the circle, and the term in the EH associated with it. 
In the limiting case when the entangling region spans 
the circle, this ’’local boundary entropy” reproduces the 
boundary entropy. 

Concretely, we consider the EH on a region V = 
U(ajbj) as depicted in Fig. [l] Consider a chiral Dirac 
fermion in on a ring with periodicity x ~ x + 2ttR. As¬ 
suming generic BC, Ik (a; + 2irR) = '&(x)e 2 '* m , the mode 

expansion is Efc e z+a b k e~ l R {t ~ x \ where 

a £ [0,1). Here a = \ for the NS sector and a = 0 
for the R sector. For a / 0, evaluating the equal time 
Green’s function on the ground state |fl), defined by 

b k > ol^) = &I< 0 |fi> = § ives 

G a (x,y) = (n\'^(x)^{y)\n)=e M r 1 n(x,y). (5) 


where n(x, y ) = (x\n\y) is the kernel of the single particle 

i ikx 

projector n onto momentum modes (x\k) = ^ e~R~ 

with k a non-negative integer |29j . Explicitly: 


n(x,y) 


1 

2nR 


E ik(x — y+iO+ ) 
e * 

k =0 


( 6 ) 


We note that the G a for a / 0 are related by spectral 
flow, implemented by a similarity transformation 

G a = U a nU~\ q/0, (7) 


where U a are unitary transformations which translate 
momenta by a: U a \k) = \k + a). Moreover, U a com¬ 
mutes with spatial projection Py, which implies that the 
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projected resolvents Q and EH kernel ([3]) for different 
a are also related by a similarity transform. Thus, we 
can relate all the a/0 resolvents L a = U a NU~ 1 to the 
following one: 


N(P) = (PynPv + P ~ §) ( 8 ) 

As an operator acting on functions defined on V, N(/3) 
is the resolvent for PynPy, which is an integrable opera¬ 
tor in the sense of USE Egg. In general, such operators are 
defined on a curve C on a manifold conformally related 
to the Riemann Sphere and have a kernel of the form 


where 


Z ± (x) = In 


IK 


i{x±ie) 

e r 


i(x±ie) 

e r 



)• 


(14) 


This form is the finite radius generalization of the result 
obtained in IT] for fermions on an infinite line. 

The EH in the a / 0 sector is obtained by applying 
([T]) to compute the resolvent and then integrating over (3 
as in ([3]) . Adapting the procedure described in jTj to our 
case, we find: 

H^f 0 = 4t T 2 e la(J1 R Ml npv(x, y)5{Z + (x ) - Z + (y)), (15) 


K{x,y) = '^f i {x)n{x,y)g i (y), (9) 

i 


where n/xyy) is a projector onto analytic functions ” in¬ 
side” E, and /?:,<??: are a set of given functions of E. It 
is well known that the resolvent, (1 + AT) -1 , is related to 
the solution of a matrix RH problem [20 . In our case 
C = S' 1 on the compactified complex plane Fig. [lj Since 
we only have a single term in A', the matrix RHP reduces 
to the scalar one and can be stated as follows. Define the 
function X (x) on S' 1 in terms of the position space ele¬ 
ments / and g by X(x) = 1 + 1/(0 — \)f(x)g(x ). The 
scalar RHP amounts to finding functions X + /X_ that 
are analytic inside/outside the unit disc satisfying 

X(x) = XZ 1 {x)X + {x), x€S\ (10) 

Equivalently, we can characterize X + / X_ as containing 
only k positive /negative momentum modes respectively. 
This leads to the crucial identities nX^n = XZ^n and 
nXZ 1 n = nXZ 1 , and to the resolvent formula p50] : 


NU 3) = 


0~\ 


1 - j-jfX^nX.g 


( 11 ) 


For our application / = g = Qy where Qy is the char¬ 
acteristic function of V (Qy(x) = 1 if x € V and is 0 
otherwise). Eq. (10) written as lnA'(x) = ln(X+(a:)) — 

ln(X_(x)) = ln(^ 3 t-) 0 y(a:), has the standard solution 
P 2 


In X± (x) = ih(0) In 


^{xyie) _ 


e* 


e ^(x=pie) _ e ibj 


( 12 ) 


where we introduced the shorthand notation 


= A In' 3 


1 

_2 
1 • 


27T 0-1 


(13) 


Substituting X± from (12) into ( flT|) , we find for the re¬ 
solvent: 


(x\N(P)\y) 


jzrj2 6 (x-y) - 


B -ih(Z+(x)-Z+( y)) 


n(x,y) 


where npy(x,y ) = n(x,y) — \d(x — y) is the principal 


part of n. Using the distribution kernel (15) to compute 


explicitly the EH, gives, owing to the S function, a local 
contribution 


^fo 0 c. = - 27r / dx 
Jy 


i d 


(1 —2a) 




2R\Z'\ V2| Z‘ 

as well as a bi-local contribution with kernel: 

(x — y ) 


*(16) 


H\ 


a^O 
V, bi-loc. 


,= 2 ^ 22 


\Z\x)\ 


-1 




R( 1-e^) 


-6{x-yi(x)), (17) 


where yi{x) are solutions of Z(x) = Z(y). 

Let us pause to inspect the single interval result for 
a = 1. Since the free fermion is a conformal field the¬ 
ory, the CFT results of uni Hu suggest that the single 
interval EH should be an integral of the energy density 
Q oc At first glance, © seems at odds with 

this form, as it contains an additional piece proportional 
to the number density p(x) = iIH(x)*(x). On closer in¬ 
spection, we see that the CFT expectation and the di¬ 
rect calculation are, in fact, consistent. Indeed, the in¬ 
terpretation of the CFT local energy tensor in mm 
as an operator, has to be done with care: in particu¬ 
lar, the form —if dx c(x)\f , t9\f r , is not Hermitian by it¬ 
self. However, integration by parts shows that it is made 
Hermitian with the addition of a term proportional to 
an integral over the density (i/2) f dx c'(x)p(x), which 
is exactly the extra term in Eq. (16) for a = 1/2. In 
the path integral formalism of m, this term can be de¬ 
rived by identifying the EH as the generator of confor¬ 
mal rotations fixing the end points of V, which must 
also generate a phase rotation of the Dirac fermion due 
to its non-trivial conformal spin. This symmetry gener¬ 
ator is the hermitian form of the local energy operator 
given by Too = Given this definition, 

our result (16) takes the form Hy = f v dx/3(x)Too, with 
0(x) = 27t|Z'(x)| _1 = 47ri?csc g ^// sin sin as a lo¬ 
cal entanglement temperature. Turning on a flux then 
introduces a shift of Too by a conserved charge density 
—pp(x) with chemical potential p, = . As shown in 
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El, this leads to a generalized first law of EE in which 
small excitations of the vacuum causes a change in EE 
given by 6S v = f v dx f3(x)5((T 00 - pp(x))) . 

The Ramond sector in which a = 0 requires sepa¬ 
rate consideration, since the zero mode k = 0 acts on 
a doubly degenerate ground state subspace of occu¬ 
pied/unoccupied modes. We choose the state dictated by 
the zero-temperature limit of the Fermi-Dirac distribu¬ 
tion. For a topological zero mode, Fermi-Dirac gives the 
mixed state \ (|occupied)(occupied| + |empty)(empty|). 
In this state fJTj, the Green’s function is 

G“=° = n-i|0)<0| (18) 


to the NS. The EE Sy = —Trpv In py can be expressed 
as in integral form: 


Sv =/rd/3TV[(/3- l/2)(L(/3)-L(-/3))-^] (23) 


It follows, using ([8]), that all Dirac fermions with a/0 
have the same EE (disregarding possible anomaly contri¬ 
butions coming from the UV cutoff). 


Using (231 we find the contribution to the EE from the 


change in BC, SS = Sr — S/vs : 


SS = ^ / 0 °° G?/itanh((^)(coth(fi7r) — 1). (24) 


Expanding in the ratio ^ gives: 


and our Resolvent is 

L“=°(/3) = ( PynPy ~\+P~ iPv|0)(0|Py)- 1 (19) 

where |0) is the k = 0 momentum mode (a;|0) = ^ttR ■ 
The difference relative to the case we have already 
solved is that the zero modes introduce a shift of the 
projector n by a rank one perturbation — 1 |0) (0|. We pro¬ 
ceed by treating the zero mode contribution as a Dyson 
perturbative expansion, which we subsequently re-sum 
(see supp. material), leading to 


r a—0 


(0) = 


N((3)Py\0)(0\PyN(f3) 
2-(0\P v N((3)P v \0) ' 


( 20 ) 


We note in passing that Eq. (|20| may be applied for 
computing the EH for single particle excitations of the 
vacuum and that the re-summation also works for higher 
rank perturbations. 


While the first term on r.h.s of Eq. (201 will yield 


a similar Hamiltonian as we considered up to now, the 
second term on the r.h.s. is different: it is due to the R 
sector zero mode and produces a non-local contribution 
to the EH. Explicitly 


(zl^ero-modeO 3 )!?/) 


2sinh 2 (-Kh)e ihlZiy ' , - ZM) 

h,h 1 

7ri?(l-fe R ) 


( 21 ) 


where l v = ~ a i ) ' s the total length of V. Using 

([3]), the contribution from T^ ro _ mode to the EH is: 


^zero-mode = =W / (22) 

— oo l + e R 

= Ez \ZUx)\R S ( X - Vl(x)) + p.V. 


21 . 


” sinh (t^( z (v)- z ( x ))} 


Thus, the zero mode induced part of the EH has a non¬ 
local contribution, even for a single interval. It is re¬ 
markable that for one interval the S(x — y) term exactly 
cancels the 1/2 shift in a (the ’’chemical potential” term) 
in (16), and we get that, as far as the strictly local terms 
are concerned, H vloc fa = 0) = H V \ 0 , : {a = 1/2). 

The non-local contributions to the resolvent due to the 
zero mode also changes the EE in the R sector relative 


~ (2 2 "-l)P 2 »C(2n)e _ ll _ . . 

^ 2n(2irR) 2n ~ 48P 2 5760P 4 ’ 1 ’ 


where P 2 n are Bernoulli numbers. Remarkably, using the 
replica trick, the authors of [26] have given the change 
in the EE induced by BC change for a massless Dirac 


fermion in the form of the series (251, and noted that the 


series is not convergent. Here we have obtained the exact 


expression Eq. (24) for the EE and can now identify the 
replica trick result of |26j as an asymptotic expansion 


of the convergent integral (24). Finally, we consider the 
limiting case where V is the entire circle. Plugging l v —► 
27tP into ( |24|) , we find SS = log2, as expected from the 
degeneracy |32|. 

Finally, we find that a similar story holds for the Majo- 
rana fermion. We derived a similar relation between the 
Green’s function G_yi(x,y) = (if{x)ip(y)) and EH (331: 


h# 1 = Ih(p v g m p v )- 1 ~i) 


(26) 


This relation differs from the Dirac one by the important 
factor of —1/2, while Gm{x,v) agrees with the Dirac 
case in the NS and R sector. In addition, since if(x) 2 is a 
constant for Majoranas, we can ignore the (x — y ) terms 
in the EH kernel. Thus one finds that in the NS sector, 
the Majorana EH has the local term 


Ky% c =ttJ dx-^-^%l)(x)d x ip(x) (27) 

while the bilocal kernel remains the same. In the 
R sector the situation is more subtle, as the Majorana 
zero mode 6o has no natural partner to create a com¬ 
plex fermion. To get the minimal, non trivial Hilbert 
space representation of the Clifford algebra containing 
bo, we have to introduce an additional anti-commuting 
operator. A physical realization of such a situation is 
the edge of a p + ip superconductor, in the presence 
of a (tt flux) vortex in the bulk (see, e.g. [18] ) - Such 
a vortex acts to change the BC on the edge into R 
type and binds a Majorana zero energy mode to its 
core. We can combine the Majorana operator for 
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the bulk-core zero mode with our edge zero mode bo, 
forming a complex/Dirac fermion operator ground state 
ag = bo + ib' 0 ) which switches between it’s associated 
“|occupied)” and “|empty)” ground states. Since | 
by the Majorana anti-commutation relations, the Green’s 
function is same as the Dirac one in (18) for any choice 
of ground state. Up to the — \ factor, the resulting R 
sector single particle EH is equal to the Dirac one, while 
the EE changes by a factor of ). 

To summarize, in this letter we provide a framework 
for the computation of EH in the presence of topological 
defects by combining RHP and finite rank perturbation 
theory. We present the first exact calculation of an EH 
on a compact space, with and without zero modes. We 
find that for a single interval the Hamiltonian consists of 
a local term in agreement with the CFT prediction for a 
single interval, as well as a non-local term associated with 
the zero mode and with the possibility of several entan¬ 
gling intervals. Note that our method also applies for an 
exact, but non-topological zero mode, however, any small 
perturbation will shift such a mode away from exact zero 
and the ground state will be unique as opposed to the 
mixed state we are considering. We also present an ex¬ 
pression, Eq. (24), representing the local entropic signa¬ 
ture of the presence of a boundary changing operator that 
reproduces the boundary entropies computed previously. 
Immediate possible applications for the method are ac¬ 
cess to the EH and EE of excited states, systems with 
non-chiral fermions and systems with several fermion fla¬ 
vors. While we concentrated here on the simple case of 
a scalar RHP, a full matrix problem appears immedi¬ 
ately in more involved situations, in particular, for non- 
relativistic fermions with a finite Fermi sea (231 . 
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SUPPLEMENTARY MATERIAL 

Entanglement Hamiltonian for Majorana Fermions 

Consider a chiral Majorana fermion in (1 + l)d with 
Lagrangian t + d x )y ;, where the spatial dimension 
is periodic x ~ x + 2irR. The canonical equal-time anti¬ 
commutation relation is {ip(x),tp(y)} = 6(x — y), and we 
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require t/)'(x) = i/j(x). Given the mode expansion 


ip{t,x) 


1 

V2^R 




(28) 


with k £ Z for the R sector and k £ Z + \ for the 
NS sector, the ground state satisfies b k \£l) = 0 for k > 
0. As a single particle operator, the equal-time Green’s 
function in the NS sector is essentially the same as in the 
Dirac case, i.e.: 


G NS (x,y) = (n\ip{x + i0 + )ip{y)\fl) = e ( ™ v) n(x, y),(29) 


where n(x,y) = (x\n\y) denotes the same projector as in 
the Dirac case. 

Given the reduced density matrix py on a region V, 
the single particle EH for Majorana fermions, Hy , is 
defined by: 

p v = Z^e-^v ; (30) 

Uy 1 = fy dxdyH^ t (x,y)ip(x)ip(y). (31) 

We now show that that H jV 1 is obtained from the Majo¬ 
rana Green’s function by: 

HP = l -\n{{PyGP v )~ 1 -l), (32) 


where Py is a projector onto V, with (x\Py\y) = 
Qy(x)5(x — y). Here Qy is the characteristic function 
for V £ S 1 , with Oy(x) = 1 if x £ V and Qy(x) = 0 
otherwise. We note that EE has been studied extensively 
also in terms of Majorana fermions (see e.g. { 22 ]) > but 
are not aware of a simple treatment of the EH as we 
present below. To derive Eq. (32), we take the Majo¬ 
rana operators to be defined on a discrete lattice, with 
{'i/.'i, %pj} = Sij, and compute the correlation functions 
from the trace formula for paired states in |28| adapted 
to the present notation: 


TcfAv+nheBuMi = ^det(l + e 2A e 2B ) (33) 

where A, B are antisymmetric matrices. We have: 

Gij = (ip^) = ±(6ij+d a logTre-^e a ^-^^)| a=0 ) 
= \5ij + \d a log ^det(l +e~ 2H v i e 2a (l i >0'l-b'X*l))| Q=0 


( , -2H{y l ) . 

\ 1+e v / 2 


(34) 


where the antisymmetry of Hy 1 was used in going to the 
last line. Inverting these relations gives (32). 


Green’s function in the Ramond sector 

In evaluating the Green’s function in the R sector, 
we are allowed to choose any linear combination of the 


ground states, or their mixed state. However, in deriv¬ 
ing the EH, we assume states that are uniquely described 
by the Green’s function through Wick’s theorem. Such 
states must have a vanishing expectation values for prod¬ 
ucts of an odd number of Majorana operators, and re¬ 
strict our choice. 

This is ensured by taking states that preserve fermion 
parity. Fermion Parity (—1) F acts by (—1 ) F ip(—l) F = 
—il>. For states p that commute with (—1) F , we have: 

-Tr(fnl)) = Tr(p{-l) F i/j{-l) F ) = Tr(p^) (35) 
=> Tr(p0) = 0. 

In the R sector, taking into account & 0 j we find 
1 1 00 

t> I I % ^ ik(x — y-t-iO ' ) r\ 

G ^ X ' y) =2^R { ~2 + ^ e -^ } = G (X ’ V) (36) 

k —0 


The relation between the resolvent and the Riemann 
Hilber problem solution. 

Here we check that the expression is indeed the 
inverse of I\. For this, one checks that 

fng - fX^nX.g - ^ ^ fngfX+'nX-g = 0, (37) 

where for simplicity we restricted to the case of i = 1. 
Here /, g , X± are understood to be operators which 
act multiplicatively in position space, with fg = (/? — 
1/2)(XT 1 A' + — 1). Substituting into the last term in the 
previous expression leads to f(n—X+ 1 nX_—n(XZ 1 X + — 
1 )X+ 1 nX_)g = 0. This equality can be checked using 
the identities 

X-'n = nX-y (38) 

nXZ 1 = nXZ 1 n, (39) 

which follow from the fact, that by their definitions, X + 
contains only positive ks and therefore multiplication by 
X + can only increases momentum while, similarly, mul¬ 
tiplication by X_ can only decrease momentum. 


Rank one perturbation and resummation 


Consider a matrix A and a rank one perturbation B = 
|s)(s| we can formally expand in A~ 1 B and then re-sum: 


{A + B) -1 = (1 + A~ l B)- 1 A~ l (40) 

= A- 1 Ay S )(£„y-i)y lA-yryiA- 1 

_ a -i _ 

l+(s|A- 1 |s) 


Substituting A = PynPy — 1/2 + /3 and |s) = yPy |0) 
then gives the desired resolvent 


L a=0 (P) = N{0) + 


N(/3)Py\0}(0\PyN(P) 
2 - (0\PyN(P)Py\0) ■ 


(41) 















Written explicitly using the kernels of the operators 
above the resolvent kernel is given by: 
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(x | ^zero-mode my) 


1 j v dzdz'{x\N(P)\z){z'\N{P)\y) 

2itR 2 ~ 2htJ v dzdz '( z \ N m z ') 



